Regge calculus from a new angle 
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In Regge calculus space time is usually approximated by a triangulation with flat simplices. We 
present a formulation using simplices with constant sectional curvature adjusted to the presence of 
a cosmological constant. As we will show such a formulation allows to replace the length variables 
by 3d or 4d dihedral angles as basic variables. Moreover we will introduce a first order formulation, 
which in contrast to using flat simplices, does not require any constraints. These considerations 
could be useful for the construction of quantum gravity models with a cosmological constant. 
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I. INTRODUCTION 

Regge calculus [lj is an elegant discrete formula- 
tion of general relativity, where space time is approx- 
imated by a piecewise flat (simplicial) manifold. Be- 
side applications in numerical relativity it has been 
used in quantum gravity as a starting point for a 
non-perturbative definition for path integral formu- 
lations 0, Q . Whereas originally Regge calculus has 
been based on length variables, newer developments 
suggest that for four-dimensional gravity other vari- 
ables might be preferable, for instance in order to de- 
fine path integral quantization or canonical formula- 
tions [3]. In particular, spin foam models are rather 
first order formulations and additionally use areas 
(and 3d angles) instead of length variables. Attempts 
to base Regge calculus on exclusively area variables 

0. 0, H, H, El have to face the difficulty, that the con- 
straints that ensure that a piecewise geometry can be 
uniquely defined, are very non-local and not known 
explicitly. As has been recently shown Q this can 
be circumvented by introducing additional variables, 
namely the 3d dihedral angles in the tetrahedra of the 
triangulation. 

Barrett [ll[ introduced a first order formulation for 
Regge calculus, using length and 4d dihedral angles 
as variables. In order to ensure that the 4d dihedral 
angles are consistent with the geometry defined by the 
length variables one constraint per 4-simplcx had to 
be added to the action. By combining [|| and [ll[ one 
can obtain a formulation with areas, 3d angles and 4d 
angles. 

Another development [l2| concerns an improvement 
of the Regge action for vacuum gravity with a cosmo- 
logical constant. Usually the cosmological constant is 
accommodated by adding a volume term to the ac- 
tion and by still using a piecewise flat triangulation, 

1. e. simplices which have flat geometry. In this way 
even the simplest solution corresponding to homoge- 
neously (maximally symmetric) curved space is only 
an approximate one compared to the continuum so- 



lution. The reason is that one uses flat simplices to 
approximate (homogeneously) curved space. Here the 
idea is to use simplices with homogeneous curvature, 
i.e. constant sectional curvature, instead, so that the 
discrete solution with homogeneous curvature is also 
an exact solution of the continuum. It has been show 
in [T3] that this allows a better representation of the 
(diffcomorphism) symmetries of general relativity and 
therefore could simplify quantization. Indeed in 3d 
the Tuarev-Viro model [lj] gives a partition function 
for homogeneously curved simplices [l4l |. whereas a 
similar quantization based on flat simplices and with 
a cosmological constant is not available yet. In gen- 
eral using a (positive) cosmological constant has the 
advantage that it can serve as a regulator in the path 
integral, as for instance the translation symmetries 
which for flat simplices would lead to divergencies are 
now compactified to the (4-) sphere. 

The geometry of simplices with and without (homo- 
geneous) curvature differs in one important way: For 
flat simplices any set of angles can specify at most the 
conformal geometry of the simplex (i.e. all the lengths 
modulo one factor). Moreover the dihedral angles of 
a flat simplex have to satisfy one constraints, namely 
that the determinant of the so-called Gram matrix 
vanishes. This is the constraint that is added in [ll[ to 
obtain a first order formulation. For homogeneously 
curved simplices, however, the set of dihedral angles 
specifies the full geometry of the simplex. Moreover 
they do not have to specify any constraints. 1 As we 
will see this will allow us to obtain formulations where 



1 Except for certain inequalities, that can be understood to 
replace the generalized triangle inequalities for the length 
variables. For a spherical simplex the angle Gram matrix, 
denned in the appendix [A] has to be positive definite. For a 
hyperbolic D-simplex the determinant of the angle Gram ma- 
trix has to be negative, and all principle D X D sub-matrices 
have to begositive definite and the cofactors positive, see for 
instance 
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the basic variables are only angles, that is either the 
3d dihedral angles or the 4d dihedral angles. Hence 
in a path integral the integration over length vari- 
ables can be replaced by an integration over angles, 
which opens up new ways for quantization, for in- 
stance for constructing the path integral measure. In 
spin foam models one usually integrates out the met- 
ric variable and is just left with a connection variable, 
which would correspond to using only the 4d dihedral 
angles. 

In the next section |TT] we will present a first order 
formulation. Here we will see that the formulation 
with curved simplices is far less complicated compared 
to working with flat simplices. In section IIIII we will 
introduce a formulation based on 3d dihedral angles. 
These have to be constrained, and we will see that 
these constraints are again relations between angles, 
this time 2d angles. Finally we discuss in section HVl a 
formulation with 4d dihedral angles, which have also 
to be constrained, this time the constraints are re- 
lations between 3d angles. We conclude in section 
IVl and summarize the necessary background on geo- 
metric relations in simplices, in particular relations 
between dihedral angles in the appendix [A"l 



II. A FIRST ORDER FORMULATION 

We will start to describe a first order formulation. 
We will see that in contrast to using flat simplices [llj 
we do not need to add any constraints to the action 
and to obtain the equations is extremely straightfor- 
ward. To use simplices with constant curvature for a 
first order formulation was suggested in [llj], but nei- 
ther an action nor any other details have been given. 

The Regge action 2 in 4 dimensions with simplices of 
constant curvature and cosmological constant A = 3k 
(and without boundary terms) is given by [l2| 



S[l e ] = ^2a t (l)e t {l) + 3 K ^2V (r (l) 



(1) 



Here we use lengths l e associated to edges as basic 
variables. The subindex t denotes triangles, at is the 
area of the triangle t and et := 2tt — YlcrDt ®t 1S the 
deficit angle associated to the triangle t. Moreover, 
6f is the 4d dihedral angle in the simplex a between 
the two tetrahedra that share the triangle t. The 
deficit angles specify the (corrected) curvature, which 
has distributional support on the triangles. The full 
curvature is given by these deficit angles and a homo- 
geneous contribution. The latter leads to a plus sign 



in front of the volume term in the action (fTJ) as com- 
pared to the continuum expression J ^Jg{R — 2A)rf x 
for the Einstein Hilbert action, where R is the (full) 
Ricci curvature scalar. 

The variation with respect to the edge length l e of 
the action |T]) gives the equations of motion 



tDe 



\ - da t 



(2) 



The variation of the deficit angles and the volume can- 
cel out due to the Schlafli identity 



3n5V a = atS9° 

tCo- 



(3) 



which holds for any variation S of the geometry of the 
4-simplex a. 

In a 3 + 1 formulation the 4d dihedral angles Of 
would specify the extrinsic curvature and can there- 
fore be taken as first order variables. Indeed this has 
been also done in [ll[ for flat simplices. As already 
mentioned in contrast to the dihedral angles in a flat 
simplex the dihedral angles in a simplex with con- 
stant curvature do not need to satisfy any constraints. 
Indeed it turns out that we do get the correct equa- 
tions of motion if we express the volume term (and the 
deficit angles) in ([TJ) in terms of the dihedral angles: 

S[Z e ,0 t CT ]=5> t (Z)e t (0)+3K^Kr(n • (4) 



Now varying with respect to l e gives the same equation 
as before J2l) 



\ - da t 



tDe 



whereas the variation with respect to 9° gives 
- at(le) + at(O a ) =0 



(5) 



(6) 



where we use the Schlafli identity to find the vari- 
ation of the volume term. Note that the equation of 
motion ([5]) accomplishes two tasks at once, firstly it 
ensures that the dihedral angles define consistently a 
simplicial geometry (that is the 4-simplices glue prop- 
erly together), and secondly that this geometry coin- 
cides with the one defined by the length variables. 3 . 



We will work with Euclidean signature. 



3 As discussed in [Till one has to take the possibility into ac- 
count that the ten areas of a 4— simplex might allow several 
length assignments. This is however only a discrete ambigu- 
ity. 
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III. A FORMULATION WITH 3D DIHEDRAL 
ANGLES 

As mentioned in the introduction spin foam models 
suggest to use areas and 3d dihedral angles instead 
of length variables as basic variables. These variables 
can be easily constructed out of gauge formulations 
for gravity, such as the Plebanski formulation [l|| on 
which spin foams are based, see also [4J . Such a formu- 
lation was developed in Q for flat simplices. The basic 
idea is to start with the geometry of the tetrahedra. 
The geometry of one tetrahedron is described by six 
lengths or equivalently by the four areas and six 3d di- 
hedral angles satisfying four Gaufi constraints. These 
constraints allow to express four of the six dihedral an- 
gles as a function of the four areas and the remaining 
two angles (which have to be non-opposite). We take 
the same starting point for the curved simplices. Here 
the difference with the flat case is that we can express 
the four areas as a function of the six 3d dihedral an- 
gles (see appendix I A So one can either decide to 
keep the areas and to add these expressions as con- 
straints or to just work with the six dihedral angles 
per tetrahedron. Following the latter route we have 
furthermore to ensure that the tetrahedra in one sim- 
plex properly glue together. That is the geometry of 
the triangle shared by two tetrahedra has to coincide 
if defined by the two sets of dihedral angles associated 
to the two tetrahedra. Again since the three angles 
in a triangle determine the geometry it is sufficient to 
ensure that these three 2d angles coincide if calculated 
from the two sets of 3d dihedral angles. Consider a 
4-simplex and label the vertices by p = 1, . . . , 5. As 
shown in the appendix the relation between 3d dihe- 
dral angles and 2d dihedral angles in a 4-simplex a 
is 



cos ai m ,kp((i>) 



COS (f>lr. 



COS (pkl,p COS <pkm,p 



sin pkLp sm<pkm, P 



(7) 



Here 4>i m , P is the dihedral angle in a tetrahedron u(p) 
opposite the vertex p and at the edge opposite the 
edge (Im). The ai„ l: k P is the 2d dihedral angle in 
the triangles a(kp) opposite the edge (kp) and at the 
vertex opposite the edge (Im). The constraints that 
have to hold between the 3d angles <fi ar e as follows: 
the angles a in the triangle a (hp) have to coincide if 
computed from either the 3d angles in the tetrahedron 
a(k) or the 3d angles in the tetrahedron <j(p), that is 



C, 



lm,kp 



iaim,kp{4>--,p) - cos Oil m ,pk(<P-,k) ■ (8) 



As we started with the geometry of the tetrahedra 
and ensured that these glue properly to 4-simplices, 
we also enforced that the 4-simplices properly glue 
together: For this the geometry of every tetrahedron 
shared by two simplices has to coincide. This is satis- 
fied by construction, as we took the geometry of the 



tetrahedra as independent variables (and constrained 
it afterwards to ensure gluing to 4-simplices). 

Finally the action for 4d Regge calculus with curved 
simplices based on 3d dihedral angles <fi an d with La- 
grange multipliers \ a e e , is 



S = 5> t (^)e t (0) + 3 K ^V CT (0) 

t CT 

E E Ke>c°A<i>) 

<y e,e'Ccr 



(9) 



where C° e , is zero, if the edges e, e' do not share a 
triangle, and coincides with Ci m ,kp for e, e' the edges 
opposite the triangles (Imk) and (Imp) respectively. 
In the appendix I A II we will give formulas for the defi- 
nition of the areas and 4d dihedral angles as a function 
of the 3d dihedral angles. Note that a priori one has to 
specify how to calculate these quantities from the 3d 
dihedral angles, e.g. for an area at one has to define 
the tetrahedron, whose set of 3d angles is used to cal- 
culate the area from. (Alternatively take the average 
over all adjacent tetrahedra.) The same holds for the 
volumes V a , that is one can first specify how to cal- 
culate lengths from the 3d dihedral angles, which can 
then in principle used to find the volume. (Again one 
possibility is to take as these lengths the average of the 
length variables calculated from the dihedral angles of 
the adjacent tetrahedra in the simplex a.) Different 
choices lead to the same result if the constraints are 
satisfied. 

The constraints ensure that one can calculate con- 
sistently length variables from the 3d angles. The in- 
verse solutions 4>(l) can be used in the action to re- 
obtain the action Hence these two action lead to 
the same equations of motions. For the same reason if 
we introduce the (first order) variables 6f and express 
the deficit angles tt and the volumes V a as functions 
of these variables, we obtain a first order formulation 
equivalent to (j4]) which uses only angles as variables. 



IV. A FORMULATION WITH 4D DIHEDRAL 
ANGLES 

Similarly one can obtain a formulation involving 
only 4d dihedral angles (and Lagrange multipliers en- 
forcing constraints). 

The ten dihedral angles of a 4-simplex determine 
uniquely its geometry and for a single 4-simplex we 
can take these 4d dihedral angles as free variables. 
If we glue the simplices together, we have to ensure 
that this can be consistently done. Two neighboring 
4-simplices might induce a priori different geometries 
for the common tetrahedron. We have to introduce 
constraints that ensure that this geometry coincides. 
To this end we use the formula expressing the 3d di- 
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hcdral angles as a function of the 4d dihedral angles 



cos t>i T 



COS Or,! cos ( 



(10) 



contributions cancel.) On the constraint hypersurfacc 
(where lengths and areas as computed from different 
simplices agree) we recover the equations of motion 



Here <pi m ,p is the dihedral angle in a tetrahedron a(p) 
opposite the vertex p and at the edge opposite the 
edge (Zm). The 4d dihedral angle 9 pm is the one at 
the triangle a(prh) opposite the edge (pm). 

The geometry of the common tetrahedron is fixed 
by the six 3d dihedral angles, hence we introduce one 
constraint per tetrahedron and edge in this tetrahe- 
dron ensuring that these 3d dihedral angles coincide if 
calculated from either of the two adjacent 4-simplices 
a and a': 



C T e (6 a , 9° ) = cos (j>l (9 a ) - cos <% (9 a ) 



(11) 



where (f>l is the 3d dihedral angle at the edge e. The 
Regge action is now 



t tot V * / a 

EE^k^') • (12) 



Here N t is the number of simplices adjacent to the 
triangle t so that the first term in (fT2|) gives the aver- 
age of the areas af attached to the same triangle but 
computed from the dihedral angles of the 4-simplices 
a adjacent to the triangle (see appendix I A Again 
the constraints allow to determine consistently the set 
of edge length so that the functions 8f(l) can be rein- 
serted into the action and one would obtain the same 
equation of motions as for the original action ([I]) . We 
can however also vary (|12p directly with respect to 
the 4d dihedral angles 9°, . Here we can again use the 
Schlafli identity and obtain 



da t 

tCa' 



) 



2n 
N't 



EE^ = o.(i3) 



Now we multiply this equation f| 1 3|) by -^pr^ i-e. the 

inverse to -^tt, where VL is the edge length of e as 

computed from the dihedral angles in a', and sum over 
H C a': 



^ di a : 

tCa' e 



T<Zcr f e£T 



dci 

oiz: 



0. (14) 



Finally summing over all a' D e! the contribution 
with the derivatives of the constraints cancel if the 
constraints are satisfied. (The sum over the simplices 
gives two contributions per tetrahedron and edge with 
opposite sign. If the constraints are satisfied these 



E da t 
dl e 

De 



e t = 



(15) 



V. CONCLUSION 

We presented different formulations for Regge cal- 
culus with cosmological constant and with simplices 
of constant sectional curvature. Using simplices with 
constant curvature leads not only to a better approx- 
imation of the continuum but leads often to more el- 
egant formulations for instance in the case of the first 
order formulation. One particular interesting feature 
is that one can take as basic variables angles, either 
the 3d dihedral angles or the 4d dihedral angles. These 
variables have to be constrained by some gluing con- 
ditions. Again curved simplices have the advantage 
that angles are sufficient to determine the geometry. 
In this way all the gluing constraints can be expressed 
by the universal relations between dihedral angles in 
n-simplices and its (n — l)-subsimplices, derived in 
appendix IA 11 

These considerations could be useful for the con- 
struction of quantum gravity models. The Tuarev- 
Viro model defines a partition function for 3d gravity 
with a cosmological constant and its semi-classical 
limit gives the Reg ge action for simplices with con- 
stant curvature [14j . In the quantum model the cos- 
mological constant is accommodated by deforming the 
SU{2) gauge group of the underlying gauge formula- 
tion to the quantum group SU(2) q . An open issue is 
whether a similar construction is possible in 4d, for 
some steps in this direction see for instance [13, [U| ■ 

To this end a gauge formulation corresponding to 
the discretized formulations presented here would be 
useful. One possibility is to consider discretizations 
of actions using Einstein-Cartan geometries [l!|, as 
the connection used there leads to the same corrected 
curvature as for simplices with constant curvature. In 
particular a formulation similar to the Plebanski ac- 
tion [16| - corresponding to using 4d and 3d dihedral 
angles and areas - could be useful to obtain spin foam 
quantizations for 4d gravity with cosmological con- 
stant. This is not only necessary to match physical 
reality but could also provide a regularization for the 
quantum gravity models (i.e. an IR cut-off for posi- 
tive cosmological constant). 
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APPENDIX A: GEOMETRIC RELATIONS IN 
SIMPLICES 

Consider a Z?-dimensional simplex in a D- 
dimensional manifold of constant sectional curvature 
k^O (i.e. the sphere S D for positive k and hyper- 
bolic space H D for negative k) consisting of D + l ver- 
tices v±,. .. , w_d+i. Denote this simplex by (123 . . . D 
1). Any subsimplex is determined by the subset 
Vi x , . . . , v in of the vertices which span this subsim- 
plex, and will therefore be denoted as (i^, ■ ■ ■ , in)- 
The subsimplices in curved space are defined to be the 
hypersurfaces with zero extrinsic curvature as embed- 
dings in the geometry of the higher dimensional sim- 
plex. These are in fact also simplices of curvature k. 
An edge (ij) is then just given by the geodesic con- 
necting Vi and vj . Denote the geodesic lengths of the 
edges (ij) by Uj. Then the (D + 1) x (D + 1) matrix 
G with entries 



where is the zj-th cofactor of G. In this way we 
can express the lengths as a function of the dihedral 
angles. For the triangle with k > we obtain 



_ cos ay 

COSVK'ij = 



cos a ik cos a, k 



sin a lk sina 7 - fe 



(A5) 



Denote the volume of the subsimplex spanned by 
+ all vertices except Vi and Vj by V^y For variations S 
of the geometry of a simplex the Schlafli identity [2l| 



(D - 1)k5V { 



12. ..D+l) 



(A6) 



holds. There is no general explicit formula (not involv- 
ing integration) for the volume of a simplex with con- 
stant curvature available for D > 3 [22j ■ For the varia- 
tion of the volume term in the action one can however 
use the Schlafli identity (|A6[) . that re-expresses this 
variation as a variation of the dihedral angles. 



G ij Cn(lij) 

where the function c K (x) is defined by 



(Al) 



c K (x) 



COS (y/~Kx) K > 

cosh (V — KX ) n < 



is called the length Gram matrix of the simplex. We 
denote by the ij-th cofactor of G, i.e. the determi- 
nant of the matrix obtained by removing the i-th row 
and j-th column of G multiplied by (— Then 
the interior dihedral angle 6*y opposite the edge (ij) 
is given by [13] 



cos 9i 



(A2) 



As the cofactor of an invertible matrix is the inverse 
multiplied by the determinant of the matrix, this for- 
mula also holds if we replace the cofactor by the in- 
verse. The length Gram matrix is invertible for non- 
degenerate simplices For a triangle, say with 
positive curvature, this reduces to 



cos a,- 



COS y/Hlij — COS \f~Klik cos y/Rl 



sin \JHlih sin ^fnl 



jk 



, (A3) 



a relation that we will encounter again between dihe- 
dral angles of different dimension. 

The angle Gram matrix is defined by Gij :— 
— cos Oij for i 7^ j and Gn — 1 . This angle Gram 
matrix and the length Gram matrix (lAlj) are in a cer- 
tain sense dual to each other. Precisely we have [2(| 



G, 



(A4) 



-« V 33 



1. Relations between dihedral angles 



We want to relate the dihedral angles 0^ in a sim- 
plex a to the dihedral angles <ft mn ,p in the subsim- 
plex o(p) not containing the vertex p and opposite the 
edge (mn). To this end we consider the Gram matrix 
Hij (p) of the subsimplex a(p) whose entries coincide 
with those of G^ for i,j ^ p. It is straightforward to 
express the inverse H M (p) of (p) using the inverse 
G kl off/,,: 



H kl (p) = G 



kl 



QpkQlp 

Gpp 



(A7) 



Inserting the definition of dihedral angles (|A2|) we ob- 
tain 



cos < 



VG u G n 



COS U„i COS 



pi COS \Jpm 



(A8) 



From IIA7D we calculate 



(~tpi f<ip 

HU (P) = G " " = G ' 1 C 1 " cos2 M ■ ( A9 ) 



Hence 



Gpp 
turns into 



cos (pin 



cos 9i m + cos 9 p i cos 8 pm 



7pl 



sin ( 



pill 



The same relation holds for flat simplices 
inversion of formulas (|A10p yields 



cos 9i m = 



cos (pi 



7n,p 



^lp,m ' 



sin 0ip, m sin Cpl m ,p 



(A10) 
. The 

.(All) 



This allows to express the 4d deficit angles as a func- 
tion of the 3d dihedral angles (if one specifies the sub- 
simplex a(p) for every dihedral angle 0/ m ). 



() 



Note that we derived the relations (|A10[) for sim- 
plices of arbitrary dimension. Hence it holds also be- 
tween the 3d dihedral angles <fii m ,p in a tetrahedron 
o{f>) and the 2d dihedral angles ai m ,kp in the trian- 
gles a(kp) in a 4-simplex a: 

COs4>lra,p + COS0fei,p COS 4>km,p / a i n\ 

cosa /mifcp = ^- — - - (A12) 

sin <pki,p sin (f>km, P 

To calculate the areas as a function of the dihedral 
angles (j) we can invoke the formula 



K«t = a.12 + uiz + a 3 i 



(A13) 



that expresses the area as a function of the three 2d 
dihedral angles ay , which in turn can be expressed as 
a function of the 3d angles with (|A12[) . Again these 



(f> can be expressed as functions of the 4d dihedral 
angles 8, see (|A10|) . so that one can express a t as a 
function of the 3d angles in a tetrahedron or the 4d 
angles in a 4-simplex. Similarly the volume of a 4- 
simplex can (in principle) either be expressed as a 
function of the ten 4d angles or the ten lengths of the 
simplex. The 4d angles and the lengths can in turn 
be expressed as functions of the and consistency of 
this procedure is guaranteed by the constraints ([8]). 
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